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Abstract 



We examine the recently observed deviation of the muon g — 2 from the 
Standard Model prediction within the framework of gravity mediated SUGRA 
models with R parity invariance. Universal soft breaking (mSUGRA) mod- 
els, and models with non-universal Higgs and third generation squark/slepton 
masses at Mq are considered. All relic density constraints from stau- 
neutralino co-annihilation and large tan /3 NLO corrections for b ^ s'j decay 
are included, and we consider two possibilities for the light Higgs: nih > 114 
GeV and nih > 120 GeV. The combined ruh, 6 — > 57 and bounds give 
rise to lower bounds on tan/3 and mi/2, while the lower bound on gives 
rise to an upper bounds on mi/2- These bounds are sensitive to Aq, e.g. for 
ruh > 114 GeV, the 95% C.L. is tan/3 > 7(5) for Aq = 0(-4mi/2), and for 
ruh > 120 GeV, tan/? > 15(10). The positive sign of the deviation implies 
/X > 0, eliminating the extreme cancellations in the dark matter neutralino- 
proton detection cross section so that almost all the SUSY parameter space 
should be accessible to future planned detectors. Most of the allowed parts 
of parameter space occur in the co-annihilation region where mg is strongly 
correlated with mi/2- The lower bound on then greatly reduces the al- 
lowed parameter space. Thus using 90% C. L. bounds on we find for 
Ao = that tan /3 > 10 and for tan /3 < 40 that mi/2 = (290 - 550) GeV 
and mo = (70 — 300) GeV. Then the tri-lepton signal and other SUSY signals 
would be beyond the Tevatron Run II (except for the light Higgs), only the 
fi and h and (and for part of the parameter space) the ei will be accessible 
to a 500 GeV NLC , while the LHC would be able to see the full SUSY mass 
spectrum. 
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The remarkable accuracy with which the muon gyromagnetic ratio can be measured 
makes it an excellent probe for new physics beyond the Standard Model. The recently 
reported result of the Brookhaven E821 experiment now gives a 2. 60" deviation from the 
predicted value of the Standard Model IQ: 

a'-P - af = 43(16) x 10"^° (1) 

where = {g^ — 2)/2. Efforts were made initially to calculate a possible deviation from 
the Standard Model within the framework of global supersymmetry (SUSY) However, 
one may show that in the limit of exact global supersymmetry, a^^^^ will vanish 0, and 
thus one needs broken supersymmetry to obtain a non-zero result. The absence of a phe- 
nomenologically viable way of spontaneously breaking global supersymmetry made realistic 
predictions for these models difficult. In contrast, spontaneous breaking of supersymmetry 
in supergravity (SUGRA) is easy to achieve, and the advent of supergravity grand unified 
models @] led to the first calculations of qSUGRA of which ^ was the first complete 

analysis. Since that time there have been a number of papers updating that result. (See 

e.g. m 

In SUGRA models, the spontaneous breaking of supersymmetry triggers the Higgs VEV 
and hence the breaking of SU{2) x f/(l), relating then these two mass scales. Thus the 
scale of the new SUSY masses is predicted to be ~ 100 GeV - 1 TeV. It was then possible 
to predict in that the SUGRA contributions would be comparable or larger than the 
electroweak contribution, 15.2(4) x 10"^*^ in accord with the now observed deviation of 
Eq. (1). This scale for the SUSY masses was further confirmed by the LEP data showing 
that consistency with grand unification could be obtained if the SUSY masses also lie in the 
above range Finally, we note that SUGRA models with R-parity invariance predict a 
dark matter candidate (the lightest neutralino) with the astronomically observed amount of 
relic density if the SUSY masses again lie in this range. 

It is thus reasonable to investigate whether the observed deviation from a^^P can be un- 
derstood within the framework of SUGRA models, and in this paper we consider gravity 
mediated SUSY breaking with R-parity invariance for models with universal soft breaking 
masses (mSUGRA) and also models with non-universal masses in the Higgs and third gen- 
eration sector. SUGRA models have a wide range of applicability including cosmological 
phenomena and accelerator physics, and constraints in one area affect predictions in other 
areas. In particular, as first observed in ||^ and emphasized in |]10|, that increases with 



tan/?, as do dark matter detection rates. Thus as we will see, the deviation of Eq. (1) 
will significantly effect the minimum neutralino -proton cross section, o"-j^o_p, for terrestial 
detectors. Even more significant is the fact that the astronomical bounds on the Xi '^elic 
density restrict the SUSY parameter space and hence the SUGRA predictions for as 
well as what may be expected to be seen at the Tevatron RUN II and the LHC. In order 
to carry out this analysis, however, it is necessary to include all the co- annihilation effects 
for large tan/?, as well as the large tan/? corrections to rrib and rrir (which are needed to 
correctly determine the corresponding Yukawa coupling constants) and the large tan (3 NLO 
corrections to the b ^ S'j decay |T^. In addition, the light Higgs {h) mass bounds play 
an important role in limiting the SUSY parameter space and it is necessary to include the 
one and two loop corrections, and the pole mass corrections. The above corrections for dark 
matter (DM) calculations were carried out in [|12|, and we will use the same corrections here. 
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Recently several papers have appeared analysing the SUGRA contribution to in light of 
the final LEP bounds on ruh and the deviation of Eq.(l) [13-16]. Relic density constraints 
were not considered in Refs. [p!5|,p!6[ and coannihilation effects apparently not included in 



Refs. [|13| , P^ . Also Refs. |[13|JT5|J16|1 do not seem to have included the constraints from the 
b —>■ s'-f decay. As will be seen below, these effects are of major importance in determining 
the SUGRA predictions. 

Before proceeding on, we state the range of parameters we assume. We take a 2cr bound 
of Eq. (1), 



11 X 10 



-10 < ^SUGRA < 75 ^ ^Q- 



10 



(2) 



a 2a bound on the b s'-f branching ratio, 1.8 x lO^'' < BR{b —>■ 57) < 4.5 x 10~^, 
and a neutralino relic density range of 0.02 < ^j^oh^ < 0.25. (Assuming a lower bound of 
0.1 does not affect results significantly.) The 6-quark mass is assumed to have the range 
4.0 GeV < mh{mh) < 4.4 GeV. We consider two bounds on the Higgs mass: rrih > 114 
GeV and rrih > 120 GeV. The first is the current LEP bound and the second is likely 
within reach of the Tevatron Run II. However, the theoretical calculations of have still 
some uncertainty as well as uncertainty in the t-quark mass, and so we will conservatively 
interpret these bounds to mean that our theoretical values obey rrih > HI GeV and 117 
GeV respectively. (Our calculations of rrih are consistent with |]l3-) The scalar and gaugino 
masses at the GUT scale obey (^0,177.1/2) < 1 TeV. We examine the range 2 < tan (3 < 40, 
and the cubic soft breaking mass is parameterized at the GUT scale by |y4o| < 4mi/2- Non- 
universal masses deviate from universality according to rnQ{l + 6) where —1 < 6 < +1. 
Other parameters are as in |jT2[ . 



We consider first the mSUGRA model, which depends on the four parameters mo, mi/2, 
Aq, tan/3 = {H2) / (Hi) ( where (if(i^2)) give rise to {d,u) quark masses) and the sign of 
the fi parameter of the Higgs mixing part of the superpotential (W = ^HiH2). The SUSY 
contribution to arises from two types of loop diagrams, i. e. those with chargino -sneutrino 
intermediate states, and those with neutralino-smuon intermediate states. The dominant 
contribution arises from the former term with the light chargino (xf)- For moderate or large 
tan P, and when (/i ± 7712)^ -C M^, one finds 



^ SUGRA 



a 



1 



m. 



tan P 

1 _ ™2 

J- . 9 



1 - 



1 + 3^ 



^2 (-^_n4-^ 



F{x) 



(3) 



where rrii = {0^1/(^0)^1/2, i = 1;2,3 are the gaugino masses at the electroweak scale and 
ac = 1/24 is the GUT scale gauge coupling constant. (One has m 



the gluino (g) mass is rrig = 
2x^(1 — x)~^lna;, where x 



m-^.) In Eq. (3), the form factor is F{x) 



m2 = 0.8mi/2: 
(1 -3x)(l -s) 



and 

-2 



The sneutrino and chargino masses being related 



to mo and mi/2 by the renormalization group equations (RGE) |T8|. (The contribution from 
the heavy chargino, reduces this result by about a third.) One finds for large mi/2 that 
F{x) = 0.6 so that decreases as l/mi/2, while for large mo, F decreases as ln(mQ)/mo 
(exhibiting the SUSY decoupling phenomena). 



Eq. (3) exhibits also the fact discussed in [|13,0] that the sign of a^^^^^ is given by 



the sign of /x. Eq. (1) thus implies that /i is positive (as pointed out in [ll4|-|l^). This then 
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has immediate consequences for dark matter detection. Thus as discussed in pO|,pT|,[T2| , for 



/i < 0, accidental cancellations can occur reducing the neutralino-proton cross section to 
below 10'^° pb over a wide range of SUSY parameters, and making halo neutralino dark 
matter unobservable for present or future planned terrestial detectors. Thus this possibility 
has now been eliminated, and future detectors (e.g. GENIUS) should be able to scan almost 
the full SUSY parameter space for mi/2 < 1 TeV. 

The lower bound of Eq. (1) plays a central role in limiting the /i > SUSY parameter 
space, particularly when combined with the bounds on the Higgs mass and the 6 — > 57 
constraints. As seen above, lowering tan/? can be compensated in by also lowering mi/2- 
However, rrih decreases with both decreasing tan/3 and decreasing mi/2- Thus the combined 
Higgs and bounds put a lower bound on tan/5. This bound is sensitive to Aq since Aq 
enters in the L-R mixing in the stop (mass)^ matrix and affects the values of the stop masses. 
We find for ruh > 111 GeV (i.e. the 114 GeV experimental bound), that tan /3 > 7 for Aq = 0, 
and tan (3 > 5 for Aq = — 4mi/2- At higher rrih the bound on tan/? is more restrictive. Thus 
for rrih > 117 GeV (corresponding to an experimental 120 GeV bound), one has tan/? > 15 
for Aq = 0, and tan/? > 10 for Aq = — 4mi/2- As the Higgs mass increases, the bound on 
tan/? increases. As discussed in p2| , p!2| , p3| , pT| , for large tan/?, the relic density constraints 
leave only co- annihilation regions possible, and these are very sensitive to the value of Aq. 
Fig. 1 exhibits the allowed regions in the itlq — mi/2 plane for tan /3 = 40, rrih > 111 GeV for 
Aq = 0, — 2mi/2, and 4mi/2 (from bottom to top). The corridors terminate at low mi/2 due 
to the 6 ^ 57 and rrih constraints. Without the constraint, the corridors would extend 
up to the end of the parameter space (mi/2 = 1 TeV). We see also that the relic density 
constraint effectively determines mo in terms of mi/2 in this region. The lower bound of 
Eq. (1), however, cuts off these curves (at the verticle lines) preventing rriQ and mi/2 from 
getting too large. Thus for large tan /?, the — 1 experiment puts a strong constraint on 
the SUSY parameter space. 

The restriction of the SUSY parameter space by the constraint affects the predicted 
dark matter detection rates. Thus the exclusion of the large mo and large mi/2 domain 
of Fig. 1 generally raises the lower bounds on the neutralino-proton cross section. In Fig. 
2 we have plotted function of m^o for tan/? = 40 for the allowed corridors for 

Aq = — 2mi/2,4mi/2 and (bottom to top). The curves terminate at high m^o due to the 
lower bound on of Eq. (1). (Note that m^o = 0.4mi/2-) Again one sees the sensitivity of 
results to the value of Aq, both for the high m^o termination point and for the magnitude of 
the cross section. Over the full range one has that o"^o_p > 6 x 10~^° pb, and hence should 
generally be accessible to future planned detectors. 
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FIG. 1. Corridors in the mg — ?n.i/2 plane allowed by the relic density constraint for tan [3 = 40, 
rrih > 111 GeV, /i > for Aq = 0, — 2mx/2i 4mi/2 from bottom to top. The curves terminate at 
low nil 12 due to the b ^ s'j constraint except for the Aq = 4mi/2 which terminates due to the mh 
constraint. The short lines through the allowed corridors represent the high mi/2 termination due 
to the lower bound on of Eq. (1). 
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FIG. 2. ^x°-p ^ function of the neutralino mass m^o for tan/3 = 40, > for 
^0 = ~2mi/2; 4mi/2) from bottom to top. The curves terminate at small m^o due to the 
h constraint for = and — 2mi/2 and due to the Higgs mass bound {mh > 111 GeV) for 

^0 = 4mx/2- The curves terminate at large ra^o due to the lower bound on of Eq. (1). 

If we reduce tan (3 , one might expect the minimum value of o"^o_p to significantly decrease. 
However, the bound then becomes more constraining, eliminating more and more of the 
high mi/2, high mo region. This is shown in Fig. 3 where the minimum value of cr^o-p is 
plotted as a function of Jtj-^^o, for tan/5 = 10, > 0, rrih > 111 GeV, for Aq = — 4mi/2 
(lower curve), Aq = (upper curve). The Aq = curve terminates at low m^o due to 
the Higgs mass bound, while the A = — 4mi/2 terminates due to the b ^ s'y constraint. 
The termination at high m^o is due to the lower bound of Eq. (1). We see that the 
parameter space is now quite restricted, and so even though tan /? is quite reduced, we find 
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> 4 X 10 pb. The co-annihilation region begins at m^o 
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FIG. 3. 

curve), Aq 

Aq = 0, and the 6 — > 57 bound for Aq 
bound on of Eq. (1). 



m^o (GeV) 

= 10, /X > 0, mh > 111 GeV for ^0 



(upper 



4:1711/2 (lower curve). The termination at low m^o is due to the bound for 

-4m]^/2- The termination at high m^o is due to the lower 



If we raise rrih and require rrih > 117 GeV ( corresponding to an experimental bound of 
120 GeV), then rrih controls the termination of the curves at low m^o. Thus for tan/5 = 40, 

215 GeVfor Ao = 0, 
4mi/2 (i- e. the Aq = 4mi/2 curve is almost completely 



the curves of Fig. 2 start at m^o = 200 GeV for Aq = — 2mi/2, at m^o 
and at m^o 



246 GeV for An 



eliminated by the rrih constraint). One has thus only a narrow range of allowed m^o . The 
allowed range becomes even narrower with decreasing tan/3, and the entire parameter space 
is eliminated when tan (3 = 10. 

We turn next to consider non-universal soft breaking models with non-universal masses 
at Mg in the third generation squarks and sleptons and in the Higgs masses: 



rrir? 



ml{l + 53) 
ml{l + (56) 



m 
m 
m 



2 

H2 
2 

tR " 

2 
II - 



--ml{l + 62y, 
ml{l + 6-r). 



(4) 



Here qi = (iL,bL) squarks, II = (z/t-,?^,) sleptons, etc. and we assume —l<Si< +1. As 
discussed in |12|, the value of /i significantly controls both the relic density and cr^o_p, and 
one may understand qualitatively how fi varies from its analytic expression which is valid 
for low and intermediate tan/5: 



/i 



t2 - 1 

l + Do 



62 + 



2 
Si 

t2 



1^ 1-Do 

H n ('^3 + <>4) 



t2 



rriQ + universal parts + loop corrections. 



(5) 



Here t = tan/5, and Dq = 1 — (mt/200 GeVsin/3)^ = 0.25. One sees that the universal rriQ 
term is quite small, and one can easily choose the 5j to make the coefficient of negative. 
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A reduction of /i^ increases the higgsino content of the neutrahno, and thus increases the 



X? 



Z couphng. In [12 



it was shown that this allowed the opening of a new region of 
allowed relic density at high mi/2 and high tan/3. We consider first the simple case where 
only ^2 is non zero and choose ^2 = 1- Fig- 4 shows cr^o_p as a function of mi/2 for this 
case when tan/3 = 40, /i > 0, m^, > 111 GeV and Aq = mi/2. The lower line corresponds to 
the usual stau-neutralino co-annihilation corridor. The upper dashed curves show the new 
allowed band arising from increased early universe annihilation through the Z s-channel 
pole. It is quite broad and has a large scattering cross section. The curves terminate at low 
mi/2 due to the b ^ sj constraint, and we have terminated the curves at the high end when 
rriQ or mi/2 exceed 1 TeV. The vertical lines are the high mi/2 endpoints due to the lower 
bound on of Eq. (1). One sees that the parameter space is significantly reduced, though 
there is still a large Z-channel band remaining. Increasing rrih increases the lower bound of 
mi/2- For rrifi > 117 GeV we find the co-annihilation (solid line) now begins at mi/2 = 510 
GeV, and the Z channel band begins at 500 GeV due to the uih constraint, leaving a sharply 
reduced region of parameter space. 



-a 



a, 



. 1 



0.01 



0.001 




400 



500 



600 



700 



800 



300 



1000 



mj (GeV) 
2 



FIG. 4. 

GeV, Ao = 



cr^o_p as a function of mi/2 ("^x? 



m 



1/2 



for (Jo 



0.4mi/2) for tan/? = 40, > 0, nih > 111 
1. The lower curve is for the fi — Xi co-annihilation channel, and the 



dashed band is for the Z s-channel annihilation allowed by non-universal soft breaking. The curves 
terminate at low mi/2 due to the 6 — > 57 constraint. The vertical lines show the termination at 
high mi/2 due to the lower bound on of Eq. (1). 

A second example of new non-universal effects is furnished by choosing 6iq { = 6^ = 64^ = 
6^) to be non-zero (as might be the case for an SU{5) or 5*0(10) model). We consider here 
5io = —0.7. In this case [1^ the fi — Xi co-annihilation corridor occurs at a much higher 
value of mo than in the universal case (i. e. for mo = 600 — 800 GeV), and is somewhat 
broadened. The Z channel band lies above it and is considerably broader. In Fig. 5 we have 
plotted function of mi/2 for the lower side of the co- annihilation corridor (lower 

curve) and for the upper side of the Z channel band (upper curve) for tan /3 = 40, /i > 0, 
Aq = mi/2 and rrih > 111 GeV. (Note that while the Z channel lies at a higher mo in the 
mo — mi/2 plane than the co-annihilation corridor, the cross section is still larger since /i^ is 
reduced.) The curves terminate at the left due to the 6 ^ 57 constraint. The vertical lines 
show the termination at high mi/2 due to the lower bound on a^, significantly shrinking the 
allowed parameter space. For rrih > 117 GeV, the Higgs mass governs the termination at 
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low mi/2, and the co-annihilation (lower curve) now begins at mi/2 = 515 GeV, and the Z 
channel (upper curve) begins at mi/2 = 520 GeV. 
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FIG. 5. '^x'i'-p ^ function of for tan/3 = 40, ;U > 0, = "7-1/2 aiid m/j > 111 GeV. 
The lower curve is for the bottom of the f\ — x? co-annihilation corridor, and the upper curve is 
for the top of the Z channel band. The termination at low 7x1^12 is due to the b ^ S'j constraint, 
and the vertical lines are the upper bound on mi^2 due to the lower bound of of Eq. (1). 

The above discussion shows that for SUGRA models, and particularly for mSUGRA, the 
data, when combined with the rrih, b ^ sj and relic density constraints have begun to 
greatly limit the SUSY parameter space. Thus the rrih and b —>■ s'y constraints determine a 
lower bound on mi/2 and hence an upper bound on a^^'^'^^, while the experimental lower 
bound on determines an upper bound on mi/2- The combined and rrih bound puts 
lower bound on tan (3 for a given value of Aq. This can be seen most clearly in Fig. 6, where 
the mSUGRA contribution to is plotted as a function of mi/2 for Aq = 0, tan/3 = 10 
(lower curve), tan/3 = 30 (middle curve)and tan/3 = 40 (upper curve). Further, most of the 
allowed mi/2 region lies in the the fi — Xi co-annihilation domain ( mi/2 ~ 350 GeV), and 
so from Fig. 1 one can see that mo is approximately determined in terms of mi/2- In Fig. 6, 
the m/j bound determines the lower limit on mi/2 for tan/3=10, while b ^ s'y determines it 
for tan/? = 40. Both are equally constraining for tan/3 =30. If we consider the 90% C. L. 
bound (a^ > 21 x 10"^°) H), one finds for Aq = that tan/3 > 10, and for tan/3 < 40 that 
mi/2 = (290 — 550) GeV, and mo = (70 — 300) GeV. This greatly constrains SUSY particle 
spectrum expected at accelerators, as can be seen in Table 1. Thus at the 90% C.L. bound 
on the tri-lepton signal will be unobservable at the Tevatron Run II since tan/? and mi/2 



are relatively large and the other SUSY particles are also beyond its reach, except for 
the light Higgs, provided m/j ~ 130 GeV [|26[. (One would need to triple the Tevatron's 
energy to see a significant part of the SUSY mass spectrum.) Only the fi and ei would 
possibly be within the reach of a 500 GeV NLC (and very marginally the xf )> while all the 
SUSY particles would be accessible to the LHC. The Brookhaven E821 experiment has a 
great deal more data that can reduce the error by a factor of about 2. When analysed, this 
would greatly narrow the predictions made here. 

One of the interesting features of Fig. 6 is that mSUGRA can no longer accommodate 
large values of a^^'"'^^. If the full E821 data should require a value significantly larger 
than 40 x 10"^", this would be a signal for the existance of non-universal soft breaking. 
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From Eq.(3) one sees that one can increase by reducing /i, and from Eq. (5) this might 
be accomphshed by non-universal soft breaking of the scalar masses (and also from non- 
universal gaugino masses at Mq-) Thus the — 2 experiment may give us significant 
insight into the nature of physics beyond the GUT scale. 

70 f , 
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2 

FIG. 6. mSUGRA contribution to function of mi/2 for vIq = 0, /i > 0, for tan/? = 10, 

30 and 40 (bottom to top). 

This work was supported in part by National Science Foundation grant No. PHY- 
0070964. 

Table 1. Allowed ranges for SUSY masses in GeV for mSUGRA assuming 90% C. L. for 
for Aq = 0. The lower value of m^^can be reduced to 240 GeV by changing Aq to -4mi/2. 
The other masses are not sensitive to Aq. 
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